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ABSTRACT 

We present results from an ongoing multiwavelength radial velocity (RV) survey of the Taurus- 
Auriga star forming region as part of our effort to identify pre-main sequence giant planet hosts. These 
1-3 Myr old T Tauri stars present significant challenges to traditional RV surveys. The presence of 
strong magnetic fields gives rise to large, cool star spots. These spots introduce significant RV jitter 
which can mimic the velocity modulation from a planet-mass companion. To distinguish between 
spot-induced and planet-induced RV modulation, we conduct observations at *~6700 A and ^2.3 /xm 
and measure the wavelength dependence (if any) in the RV amplitude. CSHELL observations of 
the known exoplanet host Gl 86 demonstrate our ability to detect not only hot Jupiters in the near 
infrared but also secular trends from more distant companions. Observations of nine very young stars 
reveal a typical reduction in RV amplitude at the longer wavelengths by a factor of ~2— 3. While we 
can not confirm the presence of planets in this sample, three targets show different periodicities in 
the two wavelength regions. This suggests different physical mechanisms underlying the optical and 
K band variability. 

Subject headings: techniques: radial velocities — planets and satellites: detection — stars: pre-main 
sequence 



1. INTRODUCTION 

The discovery of over 840 exoplanett@ in the past 
twenty years has revealed that plan etary systems ar e 
common and diverse. Pulsar planet s (jWolszcza n 1994), 
hot J upiters dMavor fc Quelozlll995[ ). and retrograde or- 
bits ([Hebrard et al.ll2011[ ) repeatedly challenge our as- 
sumptions about planet formation. After two decades 
of exoplanet discoveries, the processes underlying planet 
formation remain unclear. Lacking direct observational 
inputs, theorists must deduce formation mechanisms 
from observations of mature systems. A more preferable 
approach is to catalog the planet population around very 
young (~1 Myr old) stars while these nascent worlds still 
reside in their formation environments. 

Unfortunately, observations of pre-main sequence stars 
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are complicated by star spots, jets, accretion, and cir- 
cumstellar disks. The most successful attempts at young 
planet detection to da te have come from dire ct imaging 
surveys. For example, lLafreniere et al.l (|2010f ) presented 
evidence for a giant planet around the ^5 Myr old Up- 
per Sco solar analogue 1RXS J160929. 1-210524. They 
report the presence of an ^8 Mjtjp companion at a dis- 
tance of 330 AU based on JHK adaptive optics images 
from NIRI on Gemini North. Another d irect imaging 
detection is reported in iKraus fc Ire land (2012|). They 
present evidence for a faint point source ~16 AU from 
the 2 Myr old LkCa 15 using non-redundant aperture 
masking intcrfcrometry on Keck II. They intepret this 
source as a ~6 Mjtjp protoplanet. These observations, 
and others like them, are providing the first glimpse into 
the mechanisms of giant planet formation. 

There are, however, limitations to what can currently 
be accomplished with direct imaging. This technique can 
only detect massive planets at wide separation from their 
host stars. Furthermore, imaging does not provide a di- 
rect measure of companion mass; the derived masses are 
model-dependent. Radial velocity (RV) surveys, which 
have confirmed ~95% of the known exoplanet popula- 
tion^, obtain mass limits and are more suited to identify- 
ing planets on close orbits. In particular, "hot Jupiters" 
are readily accessible via RV observations. 

Low-mass young stars, however, are challenging tar- 
gets for traditional RV surveys. These targets are faint 
at optical wavelengths owing to late spectral types, large 
distances (>100 pc), and extinction from natal dust 
cloud s. Young stars a lso have strong magnetic fields 
(e.g., Uohns-Krulll [2007) that generate large, cool star 
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spots. Spots impact RV surveys of young stars by intro- 
ducing significant jitter w hich can mimic the RV modula- 
tion i nduced by a planet (iSaar fc Donahuel|1997l : iHatzes! 
I2002t iDesort et al.ll2007t iReiners et al.ll2010f ) . Recent at- 
tempts at detecting substellar companions in young stel- 
lar populations (10-100 Myr) via RV monitoring have 
generally been unsucc essful (e.g., iPaulson et al.l 120041 ; 
iPaulson fc Yeldal I2006T) . likely because of small sample 
sizes and the intrinsic RV variability of the targets. For 
very young, 1-3 Myr old sta rs, the RV noise proble m is 
even more acute. Recently, Ivan Evken et al.l (2012) re- 
ported observations of a possible hot Jupiter orbiting a 
2-3 Myr old weak-lined T Tauri star in Orion (P sa 0.45 
days). This interpretation is based on transit data from 
the Palomar Transit Factory (PTF) project in conjunc- 
tion with optical RV observations from the Hobby-Eberly 
Telescope (HET) and Keck I. While this is an intrigu- 
ing result, their data are dominated by spot modulation. 
Confirmation must wait for further observations, particu- 
larly in the near infrared (NIR) . The youngest confirmed 
RV planet to date is a —6 Mjup planet on a 850 day orbit 
aroun d the 100 Myr old star HD 70573 (Scti awan et al.l 
12001 . 

The success of a young star RV survey depends 
on its ability to distinguish between spot-induced and 
companion-induced RV modulation. Spot-induced RV 
variations in young stars have been identified by correla- 
tions between RVs and spectral line bisector spans (e.g., 
IHatzes et al.lll997[ ) ; planet hosts do not exhibit any such 
correlation. However, there may also be no correlation 
if the absorption lines are not spectroscopically resolved, 
i.e., the projected rotation velocity (v sin i) of the star is 
comparable to or less than the ve l ocity resolution of the 
spectrograph (IDesort et al.l 120071: IHuelamo et al.l 120081: 
IPrato et al.l l2008f K A potentially more reliable method 
for distinguishing between spots and planets leverages 
the wavelength dependence of a spot-induced RV modu- 
lation amplitude. 

The reflex motion caused by a planet affects all wave- 
lengths equally. However, the amplitude of any spot- 
induced RV v ariability will be sma ll er at longer wave- 
lengths (e.g.. IHuelamo et all 120081 : IPrato et all 120081: 
IReiners et all 2010: Ma & Gc 20J3). 

This wavelength de- 
pendence arises because of the flux-temperature scaling 
in the Rayleigh- Jeans limit of blackbody radiation; the 
contrast between a photosphere and a cooler star spot 
decreases at longer w avelengths (e.g.. iVrba et all [19861 : 
ICarpenter et al.l 120011 ). Therefore, observations in the 
visible and NIR may distinguish between stellar activity 
and a true companion by comparing the RV amplitudes 
at the two wavelengths. Furthermore, the late spectral 
types of T Tauri stars produce SEDs with peak emission 
around 1-2 /im and extinction of these partly embedded 
sources is lessened at longer wavelengths. NIR observa- 
tions can therefore lead to increased signal-to-noise and 
improved RV precision. 

NIR spectroscopy has only recently begun to play a role 
in the search for substellar companions. IMartm et all 
( 2006) combined optical and H band observations to look 
for giant planets around the young brown dwarf IP 944- 
20. They concluded that the observed optical RV mod- 
ulations were d riven by inh o moge neous surface features 
(i.e. clouds). iBlake etldl (|2007l ) used high- resolution 
K band spectroscopy to investigate the presence of gi- 



ant planets around a population of 1 dwarfs, achieving a 
precision of 300 m s" 1 . They found no evidence for com- 
panions with M sin i > 2 M.uip and P < 3 day s in their 
sample of nine targets. ISetiawan et all ()2008D reported 
the optical detection of a — 10 Mjup planet on a 3.5 day 
orbit around the 10 M yr old star TW Hydra . However, 
H band observations by Huclam o et all ([2008D revealed a 
strong wavelength dependence in the velocity amplitude. 
The low H band variability casts doubt on the presence 
of a companion and sugges ts that spots cause TW Hy- 
dra's optical RV variations. IPrato et all ([2008) observed 
DN Tau and V836 Tau — two potential planet-host candi- 
dates selected on the basis of their optical variability — in 
the K band. To within their measurement precision, they 
detected no NIR RV modulation, implying that spots 
produce the apparent optical RV variability. 

Interest in targeting M dwarfs as hosts of habitable 
planets has spurred vast improvements in NI R precision, 
bringi ng it closer to that of optical surveys. iBean et all 
(2010) report a long term precision in the H band of 
~5 m s _1 on late M dw arfs using an ammoni a gas cell 
on CRIRES at the VLT. iFigueira et ail (|2010D report a 
comparable precision, also on CRIRES, but using telluric 
absorption features as a wave length reference. B ased on 
six years of NIRSPEC data, IBlake et all (pOlOt ) report 
a precision of 50 m s _1 for their sample of K dwarfs 
and 200 m s" 1 for L dwarfs a lso through the use of tel- 
luric lines. IBailev et al.l (|2012f) also report 50 m s" 1 pre- 
cision on mid-M dwarfs, increasing to —80-170 m s _1 
for young (-10 Myr old) stars, using NIRSPEC with 
telluric-based wavelength calibration. While these re- 
sults are encouraging, these efforts are all focused on 
large (8-10 meter) aperture telescopes. Very little work 
has been done to develop precision NIR RV techniques for 
smaller aperture telescopes where more observing time 
is available to the community. We recently presented 
our technique for achieving —50 m s _1 on late-K and 
early-M d warfs using the CSHE L L NIR spectrograph on 
the IR TF ([Crockett et all 1201 ID . lAnglada-Escude et all 
(2012) also show very promising results with their suc- 
cessful efforts at obtaining 20-30 m s" 1 precision in the 
K band, also on CSHE11, using a methane isotopologue 
gas cell. 

Since 2004, we have been conducting an RV survey 
in the Taurus- Auriga star forming region in an attempt 
to catalog the giant plan et and brown dwarf population 
around very young stars (IHuerta et~l . 2008; Prato et all 
2008; Crocke ttetaTll201H IMahmud" et al. [ 120111) . In this 
paper, we present new and updated analyses on nine tar- 
gets including known spotted stars and potential young 
planet hosts; the entire sample will be discussed in a 
later paper once the entire survey is complete. In fJ2] we 
provide an overview of our observing strategy. Our data 
reduction pipeline and RV analysis techniques are pre- 
sented in SJ3J including results from RV standards used 
to quantify our long-term precision. In 21 we discuss the 
results of these analyses for each of our T Tauri targets, 
as well as updated results on the known exoplanet host 
Gl 86. Finally, in we discuss the relevance of these 
results and how they can inform future studies of the 
young planet population. 

2. OBSERVATIONS 



A Search for Giant Planet Companions to T Tauri Stars 



3 



TABLE 1 

Stellar parameters of young planet host candidates and spotted star standards 



Name 


V a 


K a 


Tcfi h 


v sin i c 


Type b 


Age b 




(mag) 


(mag) 


(K) 


(km s _1 ) 


(C/W) d 


(Myr) 


BP Tau 


12.3 


7.7 


4060 


10.9 


C 


1.8 


CI Tau 


13.0 


7.8 


4060 


10.6 


C 


2.1 


DK Tau 


12.6 


7.1 


4060 


11.5 


c 


1.0 


DN Tau 


12.5 


8.0 


3850 


8.7 


c 


0.9 


Hubble I 4 


12.0 


7.3 


4060 


12.9 


w 


0.7 


IQ Tau 


14.5 


7.8 


3780 


12.0 


c 


1.4 


IW Tau 


12.4 


8.3 


4060 


6.9 


w 


2.1 


V827 Tau 


12.2 


8.2 


4060 


19.3 


w 


2.5 


V836 Tau 


13.1 


8.6 


4060 


15.0 


c 


7.9 



a SIMBAD 

b IPalla fc Sta hlcr (2002]) 

c Glcbocki & Gnacinski (2005) 
d C = classical, W = weak-lined 

Our target list mainly consists of very yo ung (1-3 Myr) 
T Ta uri stars in the nearby (—140 pc, iKenvon et al.l 
2008) Taurus- Auriga star forming region, many of which 
we have been ob serving for several yea rs. These targets 
are selected from lHerbig fc Bell ()1988[ ) based on the fol- 
lowing criteria: (1) V < 15, (2) usini < 20 km s _1 to 
ensure relatively narrow, deep photospheric lines for ac- 
curate RV analysis, and (3) GKM stars to maximize S/N 
in the red part of the spectrum for our infrared obser- 
vations. Pre-main sequence stars of late K and early M 
spectral types will evolve to become late G, early K solar 
analogues. They offer a good compromise between hav- 
ing deep ph otosp heric CO lines needed for NIR RV mea- 
surements ( §3.3j) and being massive enoug h to offer a rea- 
sonab le chance of hosting a giant planet (pohnson et al.l 
I2010D . Of the 129 T Tauri stars, -35% are classical (ac- 
tively accreting) and —65% are weak-lined (little to no 
accretion). We also include 14 targets from the older 
(—100 Myr) Pleiades open cluster. About 20% of our 
sample are known to be close (> 0".05) binaries which 
are unresolved in our observations. Although these tar- 
gets are RV variables, the timescale is on the order of 
tens to hundreds of years and thus does not impact our 
search for short-period variability. 

We conduct observations in both the optical (—6200- 
7100 A) and the K band (—2.3 /im) to identify any 
wavelength dependence on RV variability. Our observ- 
ing strategy starts with obtaining optical observations 
of a subset of targets every night for roughly a week 
at a time. A week-long observing window matches the 
typical rotation period of T Tauri stars and falls within 
the range of companion orbital periods we are capable 
of detecting. On the basis of the optical data, we then 
group our observed targets into two categories: (1) stars 
that exhibit correlations between RV variations and line 
bisector changes, and (2) stars that show no such corre- 
lations. The RV variation s of stars in Group 1 are likely 
attributable to star spots (jHuerta et a l. 2008) while stars 
in Group 2 may harbor substellar companions. How- 
ever, since a lack of RV-bisector correlation is not suffi- 
cient for definite identification of a companion, we test 
the planet hypothesis by conducting follow-up observa- 
tions of Group 2 st ars in the NIR (|Prato et al.l 120081 : 
IMahmud et 10127)11 . 

With this strategy, we only target stars that show no 



correlation between RV and bisector span for further ob- 
servations in the NIR. The question then arises as to how 
many RV-bisector pairs are needed to identify a correla- 
tion between these two quantities for a typical spotted 
T Tauri star. To address this question, we examined 
randomly sampled subsets of various sizes from the RV 
and bisector data o f the spotted young star LkCa 19 
(|Huerta et al.l l2008f ) . Using Monte Carlo methods, we 
found a correlation between the RV and bisector data 
90% of the time at a false alarm probability (FAP) of 
0.05 when we sampled twelve observations. Hence, with 
twelve measurements per star, 1 in 20 legitimate exo- 
planet systems will be erroneously rejected by this test 
and about 1 in 10 spotted stars with clear correlations 
will mistakenly pass through to the next stage of obser- 
vations. While more observations per star would improve 
these numbers, twelve observations strikes a reasonable 
balance between the need to have a large sample of stars 
and good time-coverage per star. 

2.1. McDonald Observatory 

Visible light spectra were taken at the McDonald Ob- 
servatory 2.7 meter Harlan J. Smith telescope on th e 
Robert G. Tull Coude Spectrograph ([Tull et all [19951 ). 
The cross-dispersed coude covers a wavelength range of 
3986-9952 A. We isolated the orders with strongest S/N, 
6233-7109 A, for RV analysis. A 1".2 slit yielded a spec- 
tral resolving power of R = — 60,000. Integra- 
tion times were typically 1800 s (depending on condi- 
tions) and typical seeing was —2". We took Thorium- 
Argon (ThAr) lamp exposures before and after each 
target observation for wavelength calibration; typical 
RMS values for the d ispersion solution were —4 m s -1 
(jMahmud et al.|[2TjTT[ ). 

2.2. IRTF 

The majority of our NIR observations were obtained at 
the 3 meter N ASA Infrared Telescope Facility (IRTF) us - 
ing CSHELL (|Tokunaga et al.lll990l: [Greene et al.|[l993l) . 
CSHELL is a long-slit echelle spectrograph (1.08-5.5 /im) 
that uses a Circular Variable Filter (CVF) to isolate a 
single order onto a 256x256 InSb detector array. We used 
the CVF to isolate a 50 A segment of spectrum centered 
at 2.298 /im (order 25). This region contains numerous 
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TABLE 2 
Optical RV - bisector 

CORRELATIONS 



Target 


r 


BP Tau 


0.20 


CI Tau 


0.74 


DK Tau 


0.60 


DN Tau 


0.33 


Hubble 4 


0.80 


IQ Tau 


0.28 


IW Tau 


-0.16 


V827 Tau 


-0.88 


V836 Tau 


-0.30 



deep photosphcric absorption lines from the 2.293 /Ma 
CO v = 2-0 band head as well a rich set of predomi- 
nately CH4 telluric absorption features which we use as 
a wavelength reference. The 0".5 slit yielded a typical 
FWHM of 2.6 pixels (0.5 A, 6.5 km s _1 , measured from 
arc lamp spectra) corresponding to a spectral resolving 
power of R ~ 46,000. 

We obtained data on 56 nights between Feb 2008 and 
Feb 2012. These observations included several RV stan- 
dards, a known exoplanet host, and nine T Tauri planet 
host candidates (Table [T|). At the beginning of each 
night, we imaged twenty flat fields, each with a 20-second 
integration time, using a continuum lamp to illuminate 
the entire slit. We also imaged the same number of 20- 
second dark frames. Additionally, we imaged six Ar-Kr- 
Xe emission lines, changing the CVF while maintaining 
the grating position, to determine the wavelength refer- 
ence. All of our target data were obtained using 10" nod- 
ded pairs to enable subtraction of sky emission, dark cur- 
rent, and detector bias. Integration times for each nod 
were typically 600 seconds; fo r fainter targets we took 
multiple contiguous nod pairs (jCrockett et al.ll201ll) 

2.3. Keck 

We supplemented our CSHELL data with six nights 
of observations using NIRSPEC (IMcLean et al.l H998L 
1200(1 between Mar 2007 and Jan 2012. NIRSPEC is a 
NIR, vacuum-cryogenic, high-resolution, cross-dispersed, 
echelle spectrograph operating at the Nasymth focus on 
the 10 meter Keck II telescope. We used the N7 filter 
(1.839-2.630 jUm) with the echelle and cross-disperser 
angles set to 62.72° and 36.24°, respectively, to image 
orders 30 through 35 onto the 1024 x 1024 InSb de- 
tector. This provided a wavelength coverage of roughly 
2.157-2.550 fj,m. The 0.288" x 24" slit yielded a median 
FWHM of 2. 25 pixels (0.74 A, 9.6 km s" 1 , measured from 
arc lamp spectra) corresponding to a spectral resolving 
power of R k 31,000. To optimize simultaneous telluric 
and photosphere line coverage f H3.3p . we limited subse- 
quent analyses to order 33 (2.286-2.320 /mi). 

We observed RV standards, telluric standards, and 
over a dozen potential planet hosts in Taurus. At the 
beginning of each night, we imaged approximately ten 
flat fields (3 coadds x 30 sec) using a continuum lamp to 
illuminate the entire slit. We also imaged the same num- 
ber of dark frames (3 coadds x 30 sec). Additionally, we 
imaged a Ne-Ar-Xe-Kr calibration lamp to provide an 
initial rough estimate of the wavelength reference. All of 
our target data were obtained using an ABBA 10" nod 



TABLE 3 

RV STANDARDS 



Name 


N a bs &v 


(m s 1 ) 




CSHELL 


GJ 281 


48 


66 


HD 65277 


9 


36 


HD 219538 


13 


55 


HD 225261 


11 


61 


NIRSPEC 


GJ 281 


9 


30 


HD 10476 


5 


129 



sequence to enable subtraction of sky emission, dark cur- 
rent, and detector bias. Integration times were typically 
30-60 seconds with 2-3 coadds. 

3. DATA REDUCTION 

3.1. Visible light radial velocities 

All spectra were bias-subtracted, flat-fielded, and op- 
timally extracted u sing an IDL echelle reduction code 
(jHinkle et al.1 [2000) . We determined visible light RVs 
with a cross-correlation analysis of eight echelle orders, 
each covering about 100 A. The orders were chosen for 
their high S/N, lack of stellar emission lines, and lack 
of telluric absorption lines. For each target's cross- 
correlation template spectrum, we chose the highest S /N 
observation. Using the target itself as a template reduces 
additional uncertainty from any spectral type mismatch- 
ing. Since each order was cross-correlated against the 
corresponding order in the template, our RVs are mea- 
sured relative to one observation epoch. For each obser- 
vation, the RV and its associated uncertainty were taken 
to be the mean and standard deviations, respectively, of 
the RVs from all utilized echelle orders. 

To determine our long-term RV uncertainty, we ob- 
served several stars known to be stable at a few m s : 
107 Psc, HD 4628, r Ceti HP 65277, HD 80367 
and HD 88371 dNidever et all [20021: iButler et aLl H996l: 
iCumming et al.l [19991) . The standard deviation of these 
RV standards over eight years is ^140 m s — 1 . We 
adopt this value as our long-term systematic uncertainty. 
The final uncertainty for each observation is then 140 
m s " 1 added in quadratu re with the order-to-order scat- 
ter ([Mahmud et al.ll20111 ). 

3.2. Visible light bisector analysis 

We analyzed the visible light bisect ors to determine the 
origin of each target's RV variability ([Huerta et alJl2008t 
iMahmud et alj I2011D . We used the cross-correlation 
function from each order f H3.1[) to determine a bisector 
span (the inverse of the mean slope of the bisector); the 
mean across all orders provided a mean bisector span for 
each observation. We then computed the Pearson linear 
correlation coefficcnt, r, between the bisector spans and 
RVs for each target (Figure Q] and Table [5]) . Those tar- 
gets wit h no correlation, except a few "spotted standard 
stars" f i)4.2p . were selected for NIR follow up observa- 
tions on CSHELL. 
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-2.0 -1.5 -1.0 -0.5 0.0 
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Fig. 1. — Optical RV versus bisector span for Taurus targets with correlation coefficients. All data were acquired on the McDonald 2.7-m 
telescope. Strong correlations suggest that the optical RV modulation is spot-driven. 



3.3. NIR radial velocities 

All CSHELL and NIRSPEC o bservations were pro - 
cessed using the same pipeline (jCrockett et al.l 120111) . 
The data were dark-subtracted, flat-fielded, an d opti- 
mally extracted using an IDL implementation of iHornd 
(1986). We calculated the RVs using a spectral model- 
ing technique. This technique uses two high-resolution 
template spectra to model the stellar spectrum and the 
telluric features. For the stellar spectr um, we generated 
Next Gen stellar atmosphere models (jHauschildt et al.l 
1999) tailored to the T e ff, logg, and metallicity of our 
targets. We then used SYNTHMAG (IPiskunovl Il99l 
to create template spectra fro m the NextGen m odels 
alon g with input fro m atomic (jKupka et al.1 f2000h and 
CO (jGoorvitcijll99l line lists. 

The telluric absorption features in the K-band pro- 
vide an absolute wavelength and instrumental profile ref- 
erence, similar in concept to the iodine gas cell tech- 
nique used in high- precision optical RV exoplanet surveys 
([Butler et al.l ll996). Using the atmosphere as a "gas cell" 
lets us superimpose a relatively stable wavelength refer- 
ence, which follows the same optical path as the light 
from the science target, onto our spectra. This helps 
alleviate uncertainties introduced by variable slit illumi- 
nation, changing optical path lengths, etc. We modeled 



the telluri c features using the NOAO t elluric absorption 
spectrum (|Livingston fc Wa llace 199l|). 

In order to match the model to our observations, we 
applied a number of transformations to our templates: a 
velocity shift and a power law scaling factor, rotational 
and instrumental broadening, a second order continuum 
normalization, and a second order wavelength dispersion. 
The composite model was then binned down to the reso- 
lution of our CSHELL (or NIRSPEC) data. We used an 
IDL implementation of the Levenberg-Marquardt non- 
linear least squares fitting algorithrro to determine the 
values of the model parameters which best reproduced 
our observations. 

We used Monte Carlo techniques to estimate the errors 
in the model parameters. For each observation, we gen- 
erated 100 simulated observations based on the measured 
noise. We fit a model to each of the simulated observa- 
tions which provided us with 100 sets of parameters. We 
then calculated the mean and standard deviation of the 
100 results for each parameter which we took to be the 
final best-fit value and uncertainty, respectively. 

The RV for each nod position is the velocity shift of the 
stellar template relative to the telluric-based wavelength 

9 MPFIT, written by Craig Markwardt 
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TABLE 4 

GL 86 B PARAMETERS 



0.4 0.6 
Relative Phase 

Fig. 2. — RVs for k nown exoplanet Gl 8 6 b. Filled red circles are 
CORALIE RVs from IQueloz et ail l2000h . blue stars are CRIRES 
RVs from (Figucira et al. 2010), and open green squares are our 
K band CSHELL RVs from Nov 2008 to Nov 2011. The CSHELL 
RVs have been corrected for a long term linear drift (see Fig [3j . 
The black line is an orbit model fit to the combined datasets. The 
standard deviation of our residuals is 36 m s _1 . 




2009 



2010 



2011 



2012 



Year 



Fig. 3. — Gl 86 CSHELL RVs after subtracting the planet- 
induced modulation. The dashed line is a linear fit to the residuals 
of -0.22 m s —1 d -1 . A similar long term trend was presented by 
IQueloz ct al. (2000) and indicates the presence of a second, more 
distant, companion. 

solution. A final, nightly RV was determined by calculat- 
ing the average of the individual nod RVs, weighted by 
their uncertainties. The final uncertainity in the nightly 
RV was computed by taking the weighted standard de- 
viation of the nod RVs and dividing by the square root 
of the number of nods. 

To assess the long-term uncertainties in our observa- 
tions, we routinely observed stars known to have stable 
RVs (i.e. a v < 50 m s _1 ). We observed one RV standard, 
GJ 281, on every available night to ensure we had a con- 
sistent test spanning the entire project. When time al- 
lowed, we observed additional RV standards (HD 65277, 
HD 219538, HD 225261). Table H summarizes the RMS 
RVs for each standard. We use the RMS of GJ 281, 66 
m s -1 , as the canonical value of our long-term CSHELL 
stability since it is the only target observed consistently 
on all runs. 

GJ 281 was also observed on most nights with NIR- 
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SPEC, in addition to HD 10476 (Table |3J). We use the 
RMS of GJ 281, 30 m s _1 , to quantify our night-to-night 
uncertainties because, as with CSHELL, it is the most 
frequently observed standard with the longest time base- 
line. Using GJ 281 to assess our NIRSPEC uncertainties 
also provides a common benchmark for comparisons to 
our CSHELL observations. 



4. RV ANALYSIS 
4.1. Gl 86 



In iCrockett et al.l (|201lD , we demonstrated the ability 
of CSHELL to detect giant planets by confirming the 
exoplanet companio n around Gl 86 (jQueloz et alJ[2000t 
iFigueira et aTlf^oToh. Figure E l plo ts new CSHEL L RVs 
with those from IQueloz et alT and IFigueira et "ahl along 
with an orbital fit to all datasets. Our error bars are 
determined by adding the uncertainties from our RV al- 
gorithm in quadrature with the 66 m s _1 systematic un- 
certainty ( 33. 3p. Ou r best fit orbit (TablelJ) is very close 
to the IQueloz et al.1 solution. Our data show good agree- 
ment with previously published values exhibiting a ~36 
m s _1 standard deviation in the O — C residuals. 

The CSHELL RVs in Figure [2] have been corrected for 
a long term li near drift in this object, first reported in 
IQueloz et all (|2000f ). Figure [3] plots our RVs with the 
planet modulation subtracted and a linear fit to the resid- 
uals. A Monte Carlo analysis shows that our data are 
consistent with a 0.22 ± 0.04 m s" 1 d" 1 drift. This 
is shallower t han the 0.36 m s _1 d _1 drift reported by 
IQueloz et al.l Both datasets suggest another compan- 
ion with an orbital period on the order of hundreds of 
years (a > 20 AU). Our results demonstrate the ability 
of CSHELL to detect not only short-period Jupiter-mass 
planets but also secular trends from more distant com- 
panions. 

4.2. Discussion of individual Taurus targets 

4.2.1. Overview 

From the 56 stars on our target list with a complete 
set of twelve visible light observations, only thirteen have 
met our criteria for NIR follow-up observations. Eight 
of those currently have enough high S/N observations 
to make meaningful comparisons between the two wave- 
length bands (the remaining five will continue to be mon- 
itored in upcoming observing runs). In addition to these 
"planet host candidates", we observed a known spotted 
star, V827 Tau, to help us understand how spots will 
affect the long wavelength RV behavior. These targets, 
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TABLE 5 

Optical and NIR RV amplitudes, standard deviation, and ratios between wavelengths 



Target 


Optical 


Amplitudes 
NIR 


Ratio 


Standard Deviations 
Optical NIR Ratio 


BP Tau 


2.68 


± 


0.26 


1.70 


± 


0.26 


1.45 


± 


0.28 


0.64 


± 


0.03 


0.43 


± 


0.04 


1.43 


± 


0.15 


CI Tau 


2.80 


± 


0.30 


1.90 


± 


0.17 


1.41 


± 


0.20 


0.95 


± 


0.08 


0.63 


± 


0.05 


1.50 


± 


0.17 


DK Tau 


5.64 


± 


0.31 


3.46 


± 


0.24 


1.50 


± 


0.16 


1.58 


± 


0.07 


0.85 


± 


0.05 


1.80 


± 


0.13 


DN Tau 


2.12 


± 


0.26 


0.87 


± 


0.13 


2.27 


± 


0.42 


0.59 


± 


0.05 


0.27 


± 


0.03 


2.13 


± 


0.30 


Hubble I 4 


3.34 


± 


0.22 


1.51 


± 


0.14 


2.11 


± 


0.23 


1.07 


± 


0.04 


0.41 


± 


0.03 


2.54 


± 


0.19 


IQ Tau 


4.97 


± 


0.35 


2.59 


± 


0.54 


1.75 


± 


0.41 


1.19 


± 


0.07 


0.77 


± 


0.11 


1.48 


± 


0.23 


IW Tau 


1.40 


± 


0.18 


0.87 


± 


0.20 


1.44 


± 


0.37 


0.36 


± 


0.03 


0.25 


± 


0.04 


1.36 


± 


0.23 


V827 Tau 


5.40 


± 


0.86 


1.16 


± 


0.31 


4.24 


± 


1.13 


1.82 


± 


0.06 


0.37 


± 


0.09 


4.63 


± 


1.02 


V836 Tau 


2.90 


± 


0.42 


1.06 


± 


0.30 


2.45 


± 


0.74 


0.80 


± 


0.08 


0.35 


± 


0.09 


2.11 


± 


0.54 



along with relevant stellar parameters, are listed in Ta- 
ble [J The remainder of the sample will be discussed in 
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Fig. 4. — Optical and NIR RV amplitudes (top) and standard 
deviations (bottom). The dashed lines in each plot show optical- 
NIR ratios of one, two, and four. Targets with primarily planet- 
induced modulation should fall along the line A = 1 whereas those 
with larger A values have stronger spot signatures. 



another paper. 

For each of our targets, we calculated the peak-to-peak 
amplitude of the RV variation in both optical and NIR 
wavelengths, looked for periodicity at each wavelength, 
and fit Keplerian models to signficant periods in both 
bands. To estimate the amplitudes and their associated 
uncertainties, we generated 10 5 datasets for every target, 
replacing each RV with a Gaussian random number; the 
mean and standard deviation were set to the observed RV 
and its associated error bar, respectively. For each simu- 
lated dataset, we recorded the peak-to-peak amplitudes 
for the optical and NIR RVs {Ao and Aj, respectively) 
and the amplitude ratio between the two wavelengths 
{Ra}', we did the same for the RV standard deviations in 
the two spectral bands. We then fit Gaussian functions 
to histograms of all measured values. The mean of each 
Gaussian was taken to be the canonical value of that pa- 
rameter with error bars equal to the standard deviation 
(Tabled Figure Sj). 

To look for periodicity in our RV measurements, we 
used an IDL implem entation of the CLEAN algorithm 
(|Roberts et all 119871 ) . CLEAN attempts to remove the 
couplings between physical periods and their aliases by 
deconvolving the spectral window function from the dis- 
crete Fourier transform. We calculated the power spec- 
trum for each target using the RVs from the a) optical 
observations only, b) NIR observations only, and c) com- 
bined optical and NIR datasets. 

We also fit a Keplerian model to both the optical and 
NIR data for each target phased to significant periods 
found in the periodicity analysis. To help us quantify 
differences in the optical and NIR RVs, we started by 
fitting five free parameters (K, e, T , ui, and Vo) to the 
optical data. The period was fixed to the strongest peak 
in the optical power spectrum. Since we were primarily 
interested in how the RV amplitude changes at the longer 
wavelengths, we repeated the fit to the NIR data while 
allowing only the amplitude to vary. All other parame- 
ters were fixed at the best fit values for the optical RVs. 
For targets which showed strong periodicity in the NIR, 
we performed the additional step of reversing this pro- 
cedure: we fit a model to the NIR RVs with the period 
fixed to the peak of the NIR power spectrum. We then 
fit that model to the optical data, once again allowing 
only the amplitude to vary. For clarity, we refer to these 
four models as follows: optical Kepler fit to the opti- 
cal data (OKOD), optical Kepler fit to the infrared data 
(OKID), infrared Kepler fit to the infrared data (IKID), 
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Fig. 5. — Optical RVs (left) and NIR RVs (right) of BP Tau phased to 8.28 days. On the right, red circles are CSHELL data, green 
squares are NIRSPEC. The black line is a best- fit Kepler model. The large change in amplitude between the optical and NIR Kepler models 
indicates that the RV modulation is most likely spot-induced. 
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Fig. 6. — Optical RVs (left) and NIR RVs (right) of CI Tau phased to 6.47 days. On the right, red circles are CSHELL data, green 
squares are NIRSPEC. The black line is a best-fit Kepler model. The extremely poor fit to the NIR data indicates that the periodicity 
seen in the optical RVs is not evident in the NIR. The NIR and optical RVs are most likely caused by different physical mechanisms (see 



Fig.0. 

and infrared Kepler fit to the optical data (IKOD). 

4.2.2. BP Tau 

BP Tau is our most intensively observed target. Initial 
NIR observations strongly supported a giant planet hy- 
pothesis; continued observing has failed to uphold that. 
BP Tau exhibits significant optical RV variability with 
Aq = 2.7 km s _1 and a period of 8.28 days; the optical- 
NIR peak-to-peak RV amplitude ratio (Ra) is 1.45. The 
Kepler fits to the phase-folded RV curves (OKOD and 
OKID, Figure [5]) suggest a spot-driven interpretation of 
the RV modulation. The ratio of optical-NIR Kepler 
model amplitudes (R K ) is ~4.5. The NIRSPEC data, 
in particular, show only slight variation at about the 2a 



level. 

The optical RV period closely matches photometeric 
eriods report ed in the literature of 8.3 and 8.19 days 
Richter et aLI 119921: iPercv et al.l I2006D . High time- 
resolution photometric data from the Trans- Atlantic 
Exoplanet Survey ( TrES) also show a comparable pe- 
riod of 8.309 days (jXiao et al.1 l2012t ). However, peri- 
ods of 6.1, 7.6, and 7.7 days have al s o been identi- 
fied (ISimon et all 119901; IVrba et al.|[l986t IQsterloh etaH 
119961 respectively). The correlation with known photo- 
metric periods, the ratio of both peak-to-peak and model 
amplitudes, as well as the lack of any other significant pe- 
riodicity in the NIR RVs, all point to BP Tau's observed, 
short-period RV variation being spot-induced. 
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Fig. 7. — Same as Fig. [6] but phased to the strongest NIR period at 10.87 days. A Kepler model is unable to fit the optical data at 
this period (black line, left) but does provide a good fit to the NIR RVs. This, coupled with the 6.47 day fit and optical bisector analysis, 
suggest that the optical RVs are predominately spot-driven whereas the NIR RVs arise from a different mechanism, possibly a companion. 

CI Tau, P = 6.47 days 
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Fig. 8. — Optical RVs (left) and NIR RVs (right) of CI Tau after the best fit NIR Kepler model (Figure has been subtracted from 
both data sets. On the right, red circles are CSHELL data, green squares are NIRSPEC. The black line is a new Kepler model fit to the 
modified data. The optical RV modulation is much cleaner than Fig. [6] and the NIR RV residuals, as expected, show little variation. This 
supports our suggestion that the optical RV modulation is the superposition of both spot and companion induced variability. 



4.2.3. CI Tau 

CI Tau exhibits strong optical RV periodicity at 6.47 
days. The NIR data do not phase meaningfully to this 
period but do exhibit significant variability (Ra = 1.41, 
Figure HI). In the combined CSHELL and NIRSPEC 
data, we find a peak in the periodogram at 10.87 days 
which is not seen in the optical data. Figure [7] presents 
the same optical and NIR RVs phased to this other pe- 
riod. In this case, we freely fit all parameters in the NIR 
Kepler fit and then fit the optical data by fixing all pa- 
rameters except the amplitude. The optical fit is clearly 
nonsensical; there is no periodicity at 10.87 days in the 
optical RVs. 

Assuming the IKID models a true companion, we sub- 



tracted the planet modulation from both the optical and 
NIR RVs. Figure [8] illustrates the modified RVs phased 
to the 6.47 day period. The optical RVs demonstrate 
significant improvement. A Keplerian model to the new 
optical RVs also shows an excellent fit. In this scenario, 
the optical Keplerian parameters are most likely model- 
ing the spot modulation. The same model fit to the sub- 
tracted NIR RVs, allowing only the amplitude to vary, is 
flat. 

Although for CI Tau we have fewer than twelve optical 
observations, we see a strong correlation (r = 0.74) be- 
tween the optical RVs and bisector spans. Furthermore, 
we observed significant variability in the K band with 
a different periodicity than the optical RVs. These data 
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Fig. 9. — Optical RVs (left) and NIR RVs (right) of DK Tau phased to 4.83 days. On the right, red circles are CSHELL data, green 
squares are NIRSPEC. The black line is a best-fit Kepler model. Despite being the strongest period in the optical periodogram, the fit is 
poor at both wavelengths. 

DK Tau, P = 14.93 days 
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Fig. 10.— Optical RVs (left) and NIR RVs (right) of DK Tau phased to 14.93 days. On the right, red circles are CSHELL data, green 
squares are NIRSPEC. The black line is a best-fit Kepler model. The optical RVs show no periodicity with this phasing (as evidenced by 
the attempt at Kepler fit). The NIR RV fit is much better than the one seen at 4.83 days, though with considerable scatter. This may be 
indicative, like CI Tau, of different physical mechanisms driving the optical and NIR RV modulation. 



suggest that a heavily spotted star with a planet may ex- 
hibit optical variability that is entirely consistent with a 
spot-only model for the RV modulation. Large spots may 
mask planet-induced Doppler shifts. Modeling is needed 
to verify this possibilty (Crockett et al. in prep). 

4.2.4. DK Tau 

DK Tau exhibits large RV amplitudes in both wave- 
lengths: 5.6 km s _1 in the optical and 3.5 km s _1 in 
the NIR. An amplitude ratio of Ra = 1.5 and a lin- 
ear correlation coefficient of r = 0.6 are consistent with 
spot-induced RV variability. It is not clear, however, 
that a simple spot model is sufficient to explain the ob- 
served behavior. Figure [9] shows the optical and NIR RVs 



phase-folded to the strongest period in the optical power 
spectrum, P — 4.83 days. There is significant scatter 
in the optical RVs. After subtracting the OKOD, the 
standard deviation in the residuals is ~1 km s _1 . Fur- 
thermore, the NIR RVs do not phase to this period at 
all; the amplitude of the OKID is flat with ~1 km s _1 
residuals. 

The strongest period in the NIR RVs is, as for CI Tau 
(£ 14.2. 3j) . significantly different than that seen in the op- 
tical: 14.93 days. Figure [10] shows RV data from both 
wavelengths phased to this period. The IKID is a better 
fit Oo_c ~ 760 m s" 1 ) than the OKID, but the NIR 
RVs still show considerable scatter. The optical data do 
not phase to this period. 
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Fig. 11— Optical RVs (left) and CSHELL RVs (right) of DN Tau phased to 6.40 days. The black line is a best-fit Kepler model. The 
reduced amplitude of the NlR fit indicates that this RV variability is primarily spot-induced. 



DN Tau, P = 2.65 days 
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Fig. 12.— Optical RVs (left) and CSHELL RVs (right) of DN Tau phased to 2.65 days. The black line is a best-fit Kepler model. This 
period is the strongest seen in the NIR periodogram. The optical data do not phase to this period. The NIR fit hints at a different 
mechanism driving the NIR variability but further investigation is needed. 



4.2.5. DN Tau 

DN Tau was interpreted as a spotted star in 
iPrato et all (|2008f ) on the basis of eight NIR RVs ob- 
tained in February 2008. Revisiting these data with our 
improved NIR reduction pipeline hinted at a possibly 
more complex system so DN Tau was added back to the 
CSHELL target list for more intensive follow up. DN 
Tau shows little RV-bisector correlation (r = 0.33) and 
considerable scatter in the NIR (Ai w 870 m s -1 ). The 
amplitude ratio {Ra ~ 2.3), however, seems to again 
place this target firmly in the spot-induced category. Fig- 
ure [11] shows the optical and NIR data both phased to 
6.4 days (the strongest peak in the optical periodogram) 
along with the OKOD and OKID. The ratio of Keplerian 
model amplitudes (Rk ~ 2.7) further supports the spot 



hypothesis. 

However, interpretation of DN Tau is, similar to CI 
Tau ( t!4.2.3p . muddied by a second signficant period in 
the NIR data at 2.65 days (Figure IT2]) . To fit the NIR 
data, we followed the same approach taken for CI Tau. 
The IKOD confirms the optical power spectrum analysis: 
there is no periodicity at 2.65 days in the optical RVs. 
This is another target that will benefit from increased 
NIR scrutiny. 

4.2.6. Hubble I 4 

IMahmud et al.l (|2011l ) reported that the RV modula- 
tion of Hubble I 4 is consistent with a large polar spot 
on a nearly pole-on star. Here we refine that analysis 
with additional NIRSPEC RVs and improvements to the 
NIR spectrum modeling (Figure [13]) . 
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Fig. 13. — Optical RVs (left) and NIR RVs (right) of Hubble I 4 phased to 1.55 days. On the right, red circles are CSHELL data, green 
squares are NIRSPEC. The black line is a best-fit Kepler model. The sharply reduced amplitude in the NIR and consistent phasing between 
the two wavelengths is a strong argument for spot-driven RV modulation. 
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Fig. 14. — Optical RVs (left) and NIR RVs (right) of IQ Tau phased to 1.27 days. On the right, red circles are CSHELL data, green 
squares are NIRSPEC. The black line is a best-fit Kepler model. The NIR RVs show little variability outside their error bars, indicating 
that the optical RV variability is most likely spot-induced. 



We previously published a rotation period of 1.5459 
days arguing that the slightly stronger 2.81 day period — 
seen in a Scargle periodogram analysis — was most likely 
an alias of the true period. We based this argument on an 
observed reduction in the 2.81 day power in the CLEAN 
periodogram and on a p ublished SuperWAS P photomet- 
ric per iod of 1.5483 days ([Norton et al.ll2007t) . Kiao et al.l 
(2012) also report a photometeric period of 1.547 days 
seen in TrES data (private communication). The addi- 
tion of eight NIRSPEC observations have strengthened 
the case for a 1.54 59 day period. Revis iting the peri- 
odogram analysis of M ahmud et al.l (|2011| ) with this new 
data, the 2.81 day period no longer appears in the power 
spectrum. 



4.2.7. IW Tau 

IW Tau is another target that has comparable RV am- 
plitudes in optical and NIR wavelengths {Ra ~ 1-4) and 
shows evidence of different periods in the two bands. At 
the strongest optical period (2.95 days), the optical RVs 
fit reasonably well to an eccentric Kepler model whereas 
the NIR RVs do not phase to this period at all (Figure 
l7J|) . The strongest NIR period is 19.57 days (Figure IT?)). 
Phased to this period, the optical RVs show little period- 
icity while the NIR RVs, though plagued by significant 
scatter, show an adequate fit to a Keplerian model. 

Not only does Hubble I 4 support our argument that 
optical observations are not sufficient when looking for 
planets around young, active stars, it argues that relying 
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Fig. 15. — Optical RVs (left) and NIR RVs (right) of IQ Tau phased to 9.50 days. On the right, red circles are CSHELL data, green 
squares are NIRSPEC. The black line is a best-fit Kepler model. The 9.5 day period is the strongest seen in the NIR periodogram. The 
optical RVs show considerable scatter around the model, most likely caused by star spots. The NIR fit is driven mostly by the CSHELL 
RVs; the NIRSPEC RVs do not show the same degree of variability. 

IW Tau, P = 2.95 days 
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Fig. 16. — Optical RVs (left) and NIR RVs (right) of IW Tau phased to 2.95 days. On the right, red circles are CSHELL data, green 
squares are NIRSPEC. The black line is a best-fit Kepler model. The extremely poor fit on the right shows that the strongest optical 
period does not show itself in the NIR data. There may be a different physical mechanism behind the large variability seen in the two 
wavelengths. 



solely on NIR observations is also not conclusive. Hubble 
I 4 was the first target in our survey to exhibit resolved 
RV modulation in the K band. These data alone are 
consistent with the presence of a giant planet. However, 
comparison to the much larger variability in the optical 
data, with the same periodicity and phase, demonstrate 
that the RV modulation of Hubble I 4 is activity- induced. 

4.2.8. IQ Tau 

At a typical V magnitude of 14.5, IQ Tau is one of 
our faintest McDonald targets. The optical RVs show 
periodicity at 1.27 days with Ao ~ 5 km s" 1 . The NIR 
RVs do not phase to this period and have a lower peak-to- 



peak amplitude of 2.6 km s _1 (Figure [14]) • The reduction 
in RV variability at longer wavelengths is strong evidence 
that the optical RV behavior is spot-induced. The NIR 
power spectrum does have a peak at 9.5 days, shown 
in Figure [15] However, the IKID seems to be biased 
towards a slight dip from one data point at phase 0.1. 
The NIRSPEC RVs are not consistent with this drop 
and, in fact, show very little variation overall. 

During our observing runs at McDonald in early 2012, 
imaging at the 0.8-meter revealed dramatic variations in 
magnitude. Four months of intensive optical photome- 
tery from McDonald Observatory, Lowell Observatory, 
and the U.S. Naval Observatory Flagstaff Station (to be 
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Fig. 17. — Optical RVs (left) and NIR RVs (right) of IW Tau phased to 19.57 days. On the right, red circles are CSHELL data, green 
squares are NlRSPEC. The black line is a best-fit Kepler model. At this period, the optical RVs do not phase nor do they fit a Keplerian 
model. The NIR RVs exhibit much better phase coherence than that seen at 2.95 days. There is, however, large scatter seen in the NIR 
residuals. 
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Fig. 18.— Optical RVs (left) and CSHELL RVs (right) of V827 Tau - our "spotted standard" - phased to 3.62 days. The black line is a 
best-fit Kepler model. The drastic change in amplitude between the two wavelengths is the classic signature of a heavily spotted star. 



be published in a later paper) show extreme periodic (~6 
days) BVRI magnitude variations up to 2.5 mags. We 
have not yet determined any reasonable explanation for 
the photometeric variability of this system. 

4.2.9. V827 Tau 

We added V 827 Tau to our ta rget list as a "spotted 
standard star" (|Prato et al.ll2008l ). We have used it as a 
benchmark for how a known spotted star should behave 
at both wavelengths. We see strong periodicity in the 
optical RVs at 3.62 days. Optical RVs and bisector span 
show a strong correlation (Figure[l} with a linear correla- 
tion coefficient of r = —0.88, consistent with spot-driven 
RV modulation. Figuresl4land[T8lshow a clean separation 
between the optical and NIR ranges with Ra ~ 4.2. The 



Keplerian fits in Figure [18] also show a stark contrast. 
The OKOD has an amplitude of 2.28 km s _1 , whereas 
the amplitude for the OKID is only 0.15 km s _1 — a ra- 
tio of ~15! 

V827 Tau is another excellent demonstration of the 
need for follow-up of optical observations at longer wave- 
lengths: the optical data alone fit reasonably well to a 
Keplerian model whereas the NIR data show little varia- 
tion outside of the uncertainties. Without the NIR data 
(and bisector analysis) , it would be easy to conclude that 
the RV variability was planet-induced. 

4.2.10. V836 Tau 



V83 6 Tau was also originally presented in iPrato et al.l 
(2008) with the conclusion that, despite a lack of RV- 
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Fig. 19.— Optical RVs (left) and CSHELL RVs (right) of V836 Tau phased to 2.48 days. The black line is a best-fit Kepler model. There 
is little variability outside the error bars in the NIR RVs, suggesting that this target's RV modulation is spot-induced. 

V836 Tau, P = 3.55 days 
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Relative Phase 

Fig. 20.— Optical RVs (left) and CSHELL RVs (right) of V836 Tau phased to 3.55 days. The black line is a best-fit Kepler model. 
Phased to a longer period, the optical RVs lose phase coherence whereas the NIR RVs show a slightly better fit to a Kepler model. The 
amplitude is, however, at the same scale as the error bars. The reasonable NIR fit suggests that higher precision RV follow-up (i.e., with 
NIRSPEC) is warrented before the companion hypothesis can be ruled out. 



bisector correlation (r = -0.3), the RV modulation was 
spot-induced and not the result of an unseen compan- 
ion. The numbers presented here are updated based on 
additional NIR observations and improvements to the 
CSHELL data reduction pipeline. Despite significant op- 
tical RV variability (Ao — 2.9 km s _1 ), the NIR scatter 
is almost entirely within the CSHELL error bars. With 
R A « 2.45 and R K « 2.7 (P = 2.48 days, Figure [HI), 
the data support the spot hypothesis. 

We do see some periodicity in the NIR data at 3.55 
days (Figure l2Tfl) . We repeated the approach of CI Tau 
and DN Tau by fixing the optical model parameters to 
those determined by the NIR fit ( gXg) . The IKID fits 
the NIR data better than the OKID; the infrared model 
does not fit the optical RVs at all. However, the low 



level of NIR variability prevents us from saying anything 
definitive. This is yet another target that can greatly 
benefit from intensive NIR follow-up, particularly with 
NIRSPEC. 

5. DISCUSSION 

Our initial hypothesis that spot-induced RV modula- 
tion should decrease at NIR wavelengths was primar- 
ily based on multiwavelength photome tric surveys (e.g., 
ICarpenter et al.ll200lHVrba et al.lll986l ). However, there 
was a lack of long term NIR spectroscopic monitoring 
of T Tauri stars on which to rely. Our own spot sim- 
ulations (to be published in a later pap er, but see also 
iHuerta et alj|2008t iMahmud et al.ll201lD in addition to 
other published models (|Desort et ah 2007; Rcincr s~et al.l 
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[2010HMa fc"Gil20T^ . predict that NIR RV monitoring 
reduces the stellar jitter of young stars by a factor of 2-3. 
Data from several y oung stars — - TW Hya (Huclam o et al.l 
l2008h. V827 Tau (IPrato et all l2(¥)8h. and Hubble I 4 
(jMahmud et al.l201lh — confirm this prediction. Further- 
more, observations of Hubble I 4 demonstrate another 
important point: NIR observations alone are not suffi- 
cient. Even in the K band, some young stars exhibit RV 
amplitudes on the order of hundreds of meters per sec- 
ond with good phase coherence over several years. Only 
by comparing amplitudes in at least two widely separated 
bands can we hope to distinguish between stellar activity 
and true companions. 

Observations, however, have forced us to consider that 
a cool spot model, at times, may be too simplistic. BP 
Tau, for example, exhibits high levels of NIR RV variabil- 
ity (~1.7 km s _1 peak-to-peak) which does not phase to 
the optical data, nor shows any strong periodicity of its 
own, and c an not be explaine d solely by the presence of 
star spots. lEiroa et al.l (2002) present results from simul- 
taneous optical and NIR photometry of PMS stars which 
suggest alternate explanations. They divide their sample 
into two groups: (i) stars with similar variability trends 
in optical and NIR bands and (ii) stars with different 
trends. Those stars in (i) can be explained by the pres- 
ence of star spots or variable obscuration (such as from a 
dusty disk). Stars in (ii) are more challenging and not all 
show the expected decrease in amplitude variation with 
wavelength. In some cases the NIR variability is greater 
than that in the optical; other targets get brighter in one 
band while getting dimmer in another. This suggests 
that for some PMS stars, the optical and NIR variability 
have different underlying physical mechanisms. While 
the optical light is dominated by the photosphere, the 
NIR flux has some contribution from the circumstellar 
disk. They estimate that the disk contributes more than 
50% of the K band flux in their sample and thus con- 
clude that the NIR variability of stars with a NIR excess 
predominately tracks chan ges in the disk. 

In the case of BP Tau, iFolha & Emerson! (|1999f ) cal- 
culate a K band veiling of 0.8 ± 0.3 based on high res- 
olution spectroscopy near the 2.1661 fj,m Br7 line. This 
implies that roughly 44% of the K band light is non- 
photospheric. We therefore can not rule out the possibil- 
ity that our K band observations are impacted by either 
changes in the disk structure or shock s from the base 
of ma gnetospheric accretion columns. iMuzerolle et al.l 
(2003) present models which support this hypothesis. 
They show that the NIR excess in a sample of classical T 
Tauri stars (including BP Tau) is consistent with a single- 
temperature blackbody with T ~ 1400 K. The amount 
of excess also rou ghly correlate s with accretion luminos- 
ity. Furthermore. iNajita et ahl (|2007| ) rep ort variability 
in BP Tau's short wavelength NIR excess. iDutrev et al.1 
(2003) suggest that BP Tau may represent a class of 
transition object in the process of clearing out its disk 
based on ano malously low 12 CO to dust emission ratios. 
IDutrev et all s imaging of CO emission raises the pos- 
sibility that the CO photospheric lines used to measure 
the stellar RV may be contaminated by variable CO emis- 
sion caused by asymmetries in a warm, rotating disk. A 
disk or accretion component to the spectrum modeling 
may therefore be necessary to fully disentangle the var- 
ious possible sources of K band variability. In future 



work, we plan to combine extant photometric data from 
the McDonald Observatory 0.8 meter with our spectra 
to help clarify the underlying source of RV variability in 
not only BP Tau but all of our targets. 

CI Tau, DN Tau, and possibly V836 Tau also illus- 
trate the importance of multiwavelength observations in 
the case of heavily spotted stars. In these cases, the op- 
tical spot noise completely masks what may be smaller 
RV variability from an orbiting planet. CI Tau is partic- 
ularly interesting in this regard: the strong correlation 
between optical RVs and bisector spans (r = 0.7) leads 
to the conclusion that the optical modulation is entirely 
activity induced. Nevertheless, we see strong NIR vari- 
ability at a different period than that seen in the optical 
data. Furthermore, subtracting the NIR variability from 
the optical RVs reveals a significantly cleaner modula- 
tion. In this case, we have a target where the optical 
RVs are almost entirely spot-dominated while the NIR 
RV variations have a different origin. Here, our approach 
of filtering targets by their optical activity fails. A far 
more preferable approach is to start with NIR RV mon- 
itoring. Unfortunately, we are limited by practicality. 
A lack of community- accessible, efficient, high- resolution 
NIR spectrographs dictates our emphasis on optical RVs 
for target selection. 

Additionally, we tested the significance of the periods 
for all targets with Monte Carlo simulations. In general, 
we found that the targets with the greatest number of 
observations (e.g., BP Tau, Hubble I 4) have reasonably 
low (< 6%) false alarm probabilities (FAP). For most 
other targets, the low number of observations combined 
with a high degree of variability make the periods less 
certain. To more firmly establish the true periods, and 
eventually determine orbital parameters for likely com- 
panions, many more high-cadence NIR observations are 
needed. 

6. CONCLUSIONS 

Consistent monitoring of RV standards and known ex- 
oplanet hosts has demonstrated the long term viability of 
our techniques at a time when NIR RV observations are 
coming into their own. With a stability of 66 m s~ 1 over 
four years, over 20% of the extant exoplanet population is 
accessible to our observing methodology. Furthermore, 
this success has been demonstrated on a 3-meter class 
telescope using a twenty year old Cassegrain-mounted 
spectrograph. While parallel efforts on large telescopes 
(Keck, VLT) are highly encouraging, success on a small 
telescope, where there is greater community access, is 
arguably of even greater value. The large blocks of con- 
secutive nights required for a project like this are very 
difficult to acquire on 8-10 m class telescopes. This work 
supports the findings of the recent ReSTAR (Renew- 
ing Small Telescopes for Astronomical Research) report: 
small (< 6 meters) telescopes continue to produce inno- 
vative science that not only complements efforts on 8-10 
m class facilities but also explores niches for which large 
telescopes are not well-suited. 

While we have demonstrated the capabilities of our 
techniques, we have yet to confirm the presence of any 
planets in our young star sample. Ma yor et al.l (|2011l ) 
report that 65% of Sun-like stars harbor a planet of any 
mass with P < 100 days. Of the known RV-detected 
exoplanets with P < 100 days (such that their orbits are 
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well sampled in our survey), 4% have semi-amplitudes 
greater than the 200 m s _1 needed for a robust 3c de- 
tection with CSHELlO Since the T Tauri stars of today 
will evolve to become the solar analogues of tomorrow, 
we expect 3-4 planet detections out of the entire 143 
star sample (143 x 0.65 x 0.04) if the pre-main sequence 
planet population mirrors the main sequence one03 Cur- 
rently, only CI Tau is a strong contender; the periodicity 
seen in the other two candidates, DN Tau and V836 Tau, 
is so far only suggestive. 

The lack of planet detections to date, while likely 
impacted by small number statistics, may indicate a 
real difference between the main sequence and T Tauri 
planet populations. In situ formation of hot Jupiters is 
unlikely given both the unreasonable requirements for 
delivery of solids to the inner disk and accretion effi- 
ciency. Therefore, short-period massive planets most 
likely arrive at thei r final orbits via migration (e.g., 
iKlev fc Nelson! I2012D . If hot Jupiters are not present 
around T Tauri stars, this can place a lower limit of 1-3 
Myr on the migration timescale. In this scenario, the 
Jupiters have formed, most likely beyond the snowline, 
but have not yet migrated into short-period orbits. If 
hot Jupiters arrive via a migration mechanism with a 
timescale greater than several million years, a survey of 
young star associations over a range of ages would be 
beneficial. By identifying the age at which hot Jupiters 
arc detected, such a survey could place observational lim- 
its to the migration timescale. 

After seven years of observations, we have also yet to 
discover new short-period stellar or brown dwarf com- 
panions. Similar findings have been uncovered for brown 
dwarf companions to main sequence stars. While ~16% 
of solar analogues have close (P < 5 years) companions, 
less than 1% of those are brown dwarfs; the remainder 
are stellar (~11%) or giant pla net (~5%) companions 
(the "brown dwarf desert", e.g., iGrether fc Lineweaverl 
2006). Even before completion of our survey, this null 
result on short-period brown dwarf detections to date 
sugge sts that the desert is not evolutionary in nature 
(e.g., lArmitage fc Bonnelll l2002j ) but rather intrinsic to 
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the formation of brown dwarfs. 

T Tauri stars reside in complex environments and, as 
such, require careful observations and analyses to isolate 
the various origins of their observed variability. We may 
find that between star spots, changing effective tempera- 
tures, shocks, and disk obscuration, these young stars are 
too poorly suited to RV analysis. All other approaches, 
however, have just as many, if not more, drawbacks. Di- 
rect imaging is the only technique that is beginning to 
successfully sample the young star population, but it is 
limited to very massive planets on wide orbits. Only 
an RV survey can adequately sample the inner AUs — 
possibly the terrestrial planet forming region. The po- 
tential reward therefore justifies the intensive effort re- 
quired. While the discovery of any one planet is unlikely 
to discriminate between diverse planet formation models, 
the confirmation of a planet around any of our targets 
would be paradigm-shifting. It would represent the first 
solid data point of a hot Jupiter from the epoch of active 
planet formation and provide one of the first observa- 
tional limits on the formation timescale. 
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